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Abstract: Pair production of W bosons constitutes an important background to Higgs 
boson and new physics searches at the Large Hadron Colhder LHC. We have calculated 
the loop-induced gluon-fusion process gg — > W*W* leptons, including intermediate 
light and heavy quarks and allowing for arbitrary invariant masses of the W bosons. While 
formally of next-to-next-to-leading order, the gg W*W* — > leptons process is enhanced 
by the large gluon flux at the LHC and by experimental Higgs search cuts, and increases the 
next-to-leading order WW background estimate for Higgs searches by about 30%. We have 
extended our previous calculation to include the contribution from the intermediate top- 
bottom massive quark loop and the Higgs signal process. We provide updated results for 
cross sections and differential distributions and study the interference between the different 
gluon scattering contributions. We describe important analytical and numerical aspects of 
our calculation and present the public GG2WW event generator. 
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1. Introduction 

Vector-boson pair production provides an important background to Higgs boson searches 
in the pp ^ H ^ W*W* leptons channel at the Large Hadron Collider (LHC). Since for 
dileptonic W decays no Higgs mass peak can be reconstructed, this background cannot be 
estimated from measured data via sideband interpolation. Precise theoretical predictions 
for the irreducible VF-pair continuum background are hence crucial. 

The hadronic production of W pairs has been investigated extensively in the literature 
(see e.g. Ref. [1]). The next-to-leading order (NLO) QCD corrections to qq WW 
tvl'v' have been known for some time [2-6]. More recently also single-resonant contribu- 
tions have been included [7], and the NLO calculation has been matched with a parton 
shower [8] and combined with a summation of soft-gluon effects [9]. Electroweak correc- 
tions, which become important at large WW invariant masses, have been computed in 
Ref [10]. 

In this article we present the first complete calculation of the gluon-induced process 
gg W*W* — > leptons and study its importance as a background to Higgs searches 
in the pp ^ H ^ WW — £'^£~]^rp channel. The gluon-induced background process 
is mediated by quark loops and thus suppressed by two powers of ag relative to quark- 
antiquark annihilation. Although it formally enters only at next-to-next-to-leading order, 
the importance of the gluon-fusion process is enhanced by experimental Higgs search cuts. 
These cuts exploit the longitudinal boost and the spin correlations of the WW system to 
suppress 11/^-pair continuum production through quark- antiquark annihilation [11,12]. 
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The gluon-fusion contribution to on-shell VF-pair production, gg WW, has been 
computed in Refs. [13-15]. Here, we present a fully differential calculation of gluon-induced 
W-boson pair production and decay, gg W*W* £u£'u', including the top-bottom 
massive quark loop contribution and the intermediate Higgs contribution with full spin 
and decay angle correlations and allowing for arbitrary invariant masses of the W bosons.^ 
Partial results of this work have already been presented in Refs. [20,21]. In Ref. [20] we 
found that the contributions of the first and second quark generations enhance the NLO 
WW background prediction for Higgs searches by approximately 30%. 

In the following we describe details of our calculation, introduce the GG2WW program 
and present cross sections and differential distributions. We discuss the impact of the 
third-generation contribution and interference effects between massless and massive quark 
loop as well as signal and complete background contributions. 

2. Calculation 

2.1 Amplitude calculation preliminaries 

The calculation of the 1-loop amplitude for gg — ivl'v' is sufHciently complex that it is 
advantageous to organize Feynman amplitudes using form factors of tensor integrals, which 
are then evaluated numerically. This approach works well as long as the numerical repre- 
sentation of the amplitude is stable. When calculating cross sections for tree-level processes 
at NLO, 1-Ioop ampUtudcs arc interfered with trcc-lcvcl amplitudes. The cross section for 
our loop-induced process, however, is proportional to a squared 1-loop amplitude. Stan- 
dard loop amplitude representations will thus lead to more severe numerical instabilities. 
It is therefore advantageous to employ an algebraic approach to tensor reduction that max- 
imizes the number of cancellations that occur at the analytical level. To control the size 
of the analytical expressions it is necessary to split the amplitude into irreducible building 
blocks. Thus gauge cancellations and compensations of unphysical denominators in subex- 
pressions of the full amplitude are facilitated, and one can use standard algebraic programs 
like Maple and Mathematica to factorize and simplify the expressions. The calculation 
of the amplitude proceeds in the following steps: 

- translation of Feynman diagrams to amplitude expressions; 

- amplitude organization; 

- evaluation of amplitude expressions and algebraic reduction; 

- simplification of irreducible amplitude coefficients; 

- numerical amplitude evaluation. 

^Gluon-induced tree-level processes of the type gg WWqq are expected to be strongly suppressed [16, 
17] and have thus not been taken into account. We note that results for a similar process, gg —^ Z* Z* — » 4Z^, 
have been presented in Refs. [18, 19] including massive quark contributions, correlated decays and off-shell 
effects. 
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Before describing those steps in more detail we set up our notation and provide some basic 
definitions. 

We consider gluon-induced VF-pair production and decay and thus calculate the parton 
amplitude 

gipi, Ai) + g{p2, A2) + iiP5, -) + i^(p6, +) + i^'{P7, -) + i'{P8,+) ^ , 

where £ and £' are charged, approximately masslcss leptons of different flavour and all 
momenta are ingoing. This amplitude is related to the physical amplitude by crossing 
symmetry. Ai^2 specify the gluon helicities and p| 4 = 53^4 the virtualities of the off-shell 
vector bosons. The coupling of the gluons to the vector bosons is mediated through a 
quark loop. Although six external particles are involved in the process, at most 1-loop 
4-point functions occur in the calculation, because a pure QCD initial state couples to a 
pure electroweak final state. The contributing topologies are presented in Fig. 1. 




Figure 1: Examples for topologies contributing to gg iiyt'i/'. The Z-exchange triangle diagrams 
cancel when summed. Note that all external momenta are incoming. 

We neglect masses for the first two quark generations and all leptons, and set the 
CKM matrix to unity. The photon exchange graphs vanish due to Furry's theorem. The 
Z-exchange diagrams, however, contain an axial coupling and are proportional to (m^ — 
m^)(s3 — S4), when summed over up- and down- type contributions.^ We find that they also 
vanish for massless quarks as required by Furry's theorem and weak isospin invariance. If 
nT-u 7^ n^-d, this argument is no longer valid and the triangle graphs could contribute. 
Note that in the on-shell case one has S3 = 84 = M^, and they still vanish [14]. For 
arbitrary invariant masses S3 7^ S4, we find that the contributions from double- (Fig. la) 
and single-resonant (Fig. lb) diagrams with internal Z propagator cancel each other if the 

^Note that including the single-resonant diagrams in Fig. lb) is essential to maintain gauge invariance. 
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decay leptons are massless (as assumed). The only triangle graphs that contribute are thus 
the Higgs exchange diagrams (Fig. Ic) with amplitudes proportional to the qqH Yukawa 
couplings. The box diagrams do not involve these couplings and therefore form a gauge 
invariant subset. Since only double-resonant diagrams contribute, the amplitude f actor izes 
into W*W* production and subsequent decay mediated by the chiral fermion currents 

= ^t^(P8)7''(l-75)u(p7) ■ 

We note that the gauge-parameter dependent terms of the amplitude in gauge vanish 
for massless leptons due to current conservation. The W propagators thus simplify to 
Feynman-gauge form 

—i ai^^ 



iP-M^ + iMwTw ' 

since the Goldstone bosons do not couple to the massless external leptons. 

We calculate the contributing helicity amplitudes using the spinor formalism of Rcf. [22] . 
The complexity of the calculation is governed by the number of independent scales that 
occur, i.e. six in the case at hand. We choose the Mandelstam variables s = (pi +^2)^, 
t = (^2+^3)^ and u = (pi+Ps)"^, and the W virtualities, S3 and S4, which obey the relation 
s + t + u = ss + S4. To account for the third generation, we calculate with non-zero quark 
masses mt and mj,.^ 

Following Ref. [22], we use p2 {pi) as reference vector for the polarization vector £1 
(£2) and write 

+ 1 (2-IH1-) , 1 (1-|m|2-) 



V2 (2-|l+) ' '^^"^^ V2 (1-|2+) 

M^^l'' ./^ /1+l9-\ ' '^^^^P^^ 



V2 (1+I2-) ' '^^^^^ V2 (2+|l-) 
J3^ = (6-|m|5-) , J!^ = {8-\fi\7-), 

and obtain as projectors for the ++ and -| — helicity combinations: 



Mi 

(12) s 



- {p'tP2 + PiP^ -P1-P2 9^" - i e'''''^PlaP2p) , (2.1) 

M^- ^ _ (23) [31] tr- tr- [Ai>3i>27''] .^ 

1 2 ^13^ [32] 2s{ut-S3S4) ' ^ ' 

with tr^[r] = (tr[r] — tr[75r])/2 and the spinor inner products (ij) = (p^\p^), [ij] = 
{p^\pj), where \pf) is the Weyl spinor for a massless particle with momentum pi. We 
define the epsilon tensor by Aie^^^P'^ = tr[757'*7'^7^7'^]. The spinor prefactors in Eqs. (2.1) 
and (2.2) are pure phases and can be disregarded when calculating \M\^. 



^While keeping the full rub dependence in our calculation, wc note that the limit rub ^ is a very good 

approximation for LHC energies. The induced error is 0(mi/m?) ~ 0.1%. As nib serves as an IR cutoff 
many basis function coefficients vanish in this limit. The non-zero coefficients simplify considerably, too. 
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2.2 Symbolic amplitude evaluation with algebraic tensor reduction 

As discussed in Section 2.1, the amplitude iM factorizes into the production of two virtual 
charged vector bosons 

9{pi,M) + 9{P2, A2) + W-*{p3) + W+*{p4) ^ 

and tlieir decay. The production amplitude is tlius contracted with the vector boson prop- 
agators and lepton currents: 

The coupling constants and colour factors are conveniently absorbed in the scattering tensor 
_y^/iiM2M3M4^ which can be decomposed in terms of Lorentz tensor structures built from the 
metric (?^^, the external momenta p^, P2, and the vector kg = i e^PiP2P3_ Note that kg 
products are reducible: kgk^ = ag^'^ + Pijp'^Pj- The basis of tensor structures that can 
occur is defined by momentum conservation, Schouten identities, the transversality /gauge 
conditions £1 -pi = £1 ■p2 = £2 -Pi = £2 ■P2 = and current conservation J3 -p^ = J4 -p^ = 0. 
We find 

jvi — Jig g '-'JinnjiPji Pj2Pj3 Pjt 

31, 32, j 3, j 4 

~ J3J4 y ^33 ^34 ~ 3234. y i^3l ^34 ~ 3233 y i^32 ^33 

33,34 32,34 32,33 

+ E ^3.34 9'^'^ + E ^.3 9'^'' V^V^ + Y.^.32 9'^'^ Pf. Pf. 

31,34 31,33 31,33 

+ Z^^j200'^o 9 Pj2 +Z^-^jiOo'«o 9 Pji 

32 31 

+ / .-^00 74 ^0 9 Pj^ + / .-^00.73^0 9 Pj3 

34 33 

+ ^323334'^0 Pn P33 P34 + ^3l3334'^0 Pji Pj3 P34 

32,j3,j4 31,33,34 

+ E ^1323J^^1^^1+ E 4.2.3^^^5^ P^^P^ (2-3) 



J1J2J4 31,32,33 



with Ji,j2 = 3 and j3,j4 G {l5 2}. Based on Eq. (2.3) a gauge invariant representation 
with 36 coefficients can be derived. The terms involving an epsilon tensor, i.e. k^, are 
parity odd and their coefficients are proportional to (m| — ml) and vanish if weak isospin 
is conserved [13,14], in particular for massless quarks. 

The amplitude is invariant under exchange of the gluons (Bose symmetry). Since the 
amplitude is also CP invariant,^ only two helicity amplitudes are independent: 

M—^'''\s,t,u,S3,S4) = M++^'''^{s,u,t,Si,S3), (2.4) 
M-+^^''\s,t,u,S3,S4) = M+-^^''^{s,u,t,Si,S3), (2.5) 



*Note that lepton flavour cannot be distinguished in the hmit of massless leptons. 
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with A4^1^2./3./4 = ^^1^^^2 ^t^upaj^^J^^^ 

Further simphfication can be achieved by expressing the heUcity amphtudes 

9 

^AiW4 ^ ^c}^^^s,t,u,ss,S4,mlmj) r,(J3, J4) 

in terms of the nine gauge-independent scalar structures 

Tj e { J3 • J4, J3 -Pi Ja- Pk, J3 -Pi Ji- ko, J3 -ko Ji- pk} , 

where l,k G {1,2}. The coefHcients Cj^^'^ are Hnear combinations of the amphtude coeffi- 
cients defined in Eq. (2.3) and will subsequently be expressed in terms of basis integrals. 
The final coefHcients contain negative powers of the Gram determinant 

det G = 2s{tu — S3S4) , 

which emerges during tensor reduction. In reference frames with back-to-back gluons, it is 
related to the transverse momentum of the W boson: detG = 2s^p|y. As the transverse 
momentum of the vector boson approaches zero, the inverse Gram determinant diverges 
while the amplitude remains finite. In this phase space region large numerical cancellations 
occur and finite machine precision can lead to instabilities during evaluation. To mitigate 
this effect, our goal is to cancel as many powers of det G as possible. For this purpose, we 
expose the Gram determinant in our expressions by introducing the auxiliary vector 

[detG , , , ,7, 
P3 = y (0> sm 03, cos (p3, 0) . 

For example, if ps is replaced by in Eq. (2.2) the Gram determinant in the denominator 
cancels explicity. More generally, we write the helicity amplitudes 

9 _ 

^XMJ, = ^C^^^^{s,t,u,S3,S4,mlm^t) f,(J3, J4) (2.6) 
in terms of the following scalar structures (see Appendix A): 

fj G {J3 ■ J4, Pl ■ J3P1 • J4, Pi ■ J3P3 ■ J4, Pi ■ J4P3 ■ J3, P3 ■ J3P3 - J4, e(pi,P2, J3, Ja), 
Pi - J3^{Pl,P2,P3,J4), Pl ■ J4e{pi,P2,P3,J3), P3 ' J3 ^{Pl,P2,P3, J4)} ■ 

The coefficients in Eq. (2.6) involve tensor (loop momentum) integrals, which can be 
written in terms of Lorentz tensors with coefficients involving only scalar integrals. As 
an alternative to the standard methods of Refs. [23,24], we also applied the improved 
reduction formalism of Refs. [25-27] to calculate the amplitude. Here, for instance, a rank 
two 4-point tensor integral is reduced via 

ir{ri,r2,rs, mj, m^, mi, ml) = B^'^g^'^ + ^fAh ' 

J1J2 = 1 
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where ri = pi, r2 = pi + P2 and = — p4 with the external momenta Pi,P2,P3iP4- The 
form factors -B^'^ and A'^'^j^ depend on scalar integrals and are ultimately functions of the 
Lorentz invariants Sij = {pi +Pj)^, Sj = pj and the internal masses. In general, at most 
rank three tensor box integrals can occur. ^ For the calculation at hand, we generate explicit 
analytical representations in terms of the scalar integral basis^ 

he{l!, 1^,1^,1}. 

In total 27 different scalar integrals appear. As analytical results do not exist for all 6- 
dimensional 4-point functions /f, we have represented them in terms of 3- and 4-point 
functions in n = AJ 

Since the process is loop induced, no real corrections exist and even for massless inter- 
nal quarks each loop diagram is infrared (IR) finite. We confirmed that the coefficients of 
IR divergent basis functions vanish. Since no ggWW counter term exists, the amplitude 
is also ultraviolet (UV) finite. Note however that each box diagram is UV divergent and 
only the gauge invariant sum of all box graphs is finite. We therefore employ dimensional 
regularization to define the amplitudes for individual graphs and evaluate the 2-point func- 
tion coefficients to order e = (ra — 4)/2 for space-time dimension n. The V — A coupling 
of the charged vector bosons to the internal quarks requires a prescription to treat n- and 
4-dimensional objects consistently. We apply standard dimension splitting rules [28] (see 
Appendix B). 

The coefficients in Eq. (2.6) can now be written in terms of basis functions Ik as 

Cj'^-^''{s,t,u,S3,S4„ml,mt) = ^C^l^^s,t,u, ss, S4,ml,'mt)Ik , 

k,l 

where the coefficient Cj^^''^ corresponds to diagram / and is a rational polynomial that is 
computed with and saved as Form [29] code.^ The irreducible amplitude coefficients are 
then simplified with Maple. First, each C^ll"^ fj expression is simplified. For A1A2 = 

-|— |-, all inverse Gram determinants cancel in this step. For A1A2 = H — . — + • however, 

one inverse power survives. Next, we sum over the diagrams, which facilitates further 
simplification, because discrete symmetries and gauge invariance are restored. The final 
output is converted to FORTRAN code. All steps are automatized. 

For the triangle topologies with Higgs exchange we thus obtain the well-known result 

iM{gg ^ w+^W* ^vv'l') = M^^^^-^^^^ ^ 

s — rrijj + i ruH^H 

with 

^XiX2J3J4 ^ C^^^^{s,ml,ml) J3 • J4 

^Analytical and numerical representations for the form factors are provided in Ref. [25]. 

®The tadpole integral does not appear in this list, as it can be viewed as a degenerate 2-point integral. 

'^This introduces an inverse Gram determinant, but the expression can be grouped such that the combi- 
nation of scalar integrals tends to zero as the inverse Gram determinant diverges and no additional stability 
problem is introduced. 

^Unevaluated amplitude expressions for each contributing Feynman graph have been compared with 
output from FeynArts 3.2 [30]. 
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and 



C±±(s,mg,m?) = -2m2 (s - 4m2) 7J(s, 0, 0, m?, mf , mf ) 

-2 ml {s - 4ml) 0, 0, m^, mg, mg) + 4 (m? + mg) , 
C^'^is^mlml) = 0. 

Here, only the helicity combinations A1A2 = ++ and contribute, since the intermediate 

scalar forces the gluons to be in an L = S = state. Eq. (2.4) implies C++ = C~~ . The 
explicit results for the box topologies are too complex to be presented. 

Two independent calculations of the amplitude have been performed. All symmetry 
relations of the amplitude have been checked analytically. In one calculation the tensor 
reduction methods described in this section were applied, while standard methods [23, 24] 
were used in the other. 

2.3 Numerical amplitude evaluation 

The symbolic evaluation method described in Section 2.2 strongly reduces the destabiliz- 
ing effects of inverse Gram determinants in the final amplitude representation, but does 
not completely remove them for the H — and — + helicity combinations. When evaluated 
in double precision, our analytic expression for the amplitude thus exhibits numerical in- 
stabilities in the extreme forward scattering region, e.g. pt{W^) — 0, where (detG)~^ 
diverges. Since the u pair is not detected, this phase space region still contributes to the 
cross section after application of the selection cuts. 

Numerical instabilities can be remedied by evaluating the amplitude in quadruple 
precision. But, a huge runtime penalty is incurred in comparison to double precision. In 
order to overcome this practical problem, one can restrict the use of quadruple precision 
to a small region in phase space where 

PxiW^) < 6 GeV 

and 



Pr(W+) < 1 GeV or max 



Pw± - Mw 



while double precision is used in the remainder of the phase space. Using this "mixed" 
mode of our numerical program, the results presented below in Section 3 were calculated 
with no indication of numerical instabilities and no significant runtime overhead. For a 
specific phase space configuration we compared numerical results for |A^p obtained with 
our independent amplitude calculations and found agreement. We use LoopTools [31] to 
evaluate the scalar integrals numerically. 

2.4 Cross-section calculation and the GG2WW program 

The cross sections and distributions presented in Section 3 were verified with two indepen- 
dent phase space and Monte Carlo integration implementations. 

Our public program, named GG2WW, includes all background and signal contributions, 
full spin correlations, off-shell and interference effects, as well as finite top and bottom quark 
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mass effects. It can be used either at tfie parton level or to generate weighted or unweighted 
events in Les Houches Accord format [32]. A combination of the multi-channel [33,34] 
and phase space-decomposition [35,36] Monte Carlo integration techniques was used with 
appropriate mappings to compensate peaks in the amplitude. In addition, automatized 
VEGAS-style [37] adaptive sampling is employed using OmniComp-Dvegas, which features 
a parallel mode (including histogram filling) [38] . Parton distribution functions are included 
via the LHAPDF package [39]. Selection cuts and histograms can be specified in a user- 
friendly format. The program is available on the Web [40] and has already been used by 
ATLAS and CMS in recent H WW studies [41-43]. 

3. Results 

In this section we present numerical results for the process pp — > W*W* Ivl'v' at 
the LHC. We tabulate the total cross section and the cross section for two sets of ex- 
perimental cuts. We focus on the impact of the massive top-bottom loop, which has 
been neglected in Ref. [20], and the size of the signal-background interference. As dis- 
cussed in Section 2.4, we provide a public parton-lcvel event generator for the process 
gg W*W* leptons [40], which can be used to study alternative sets of cuts or to 
generate any kind of distribution. 

The experimental cuts include a set of "standard cuts" [6], where we require both 
charged leptons to be produced at pT,e > 20 GeV and \r)e\ < 2.5 motivated by detector 
coverage, and a missing transverse momentum > 25 GeV characteristic for Icptonic W 
decays. Cross sections calculated with this set of cuts will be labeled crgfii. Various further 
cuts have been proposed for the experimental Higgs searches to enhance the signal-to- 
background ratio [11,12,44-47]. As in our previous publication [20] we have studied a set 
of cuts similar to those advocated in a recent experimental study [47]. In addition to the 
"standard cuts" defined above, we require that the opening angle between the two charged 
leptons in the plane transverse to the beam direction should satisfy A(j)T/e < 45° and that 
the dilcpton invariant mass be less than 35 GeV. Furthermore, the larger and smaller 
of the charged lepton transverse momenta are restricted as follows: 25 GeV < PT,mm and 
35 GeV < pT,max < 50 GeV. Finally, a jet veto is imposed that removes events with jets 
where PT,jet > 20 GeV and \r]jet\ < 3. Cross sections evaluated with the Higgs selection 
cuts will be labeled abkg. 

To obtain numerical results we have used the same set of input parameters as in Ref. [20] : 

Mw = 80.419 GeV, Mz = 91.188 GeV, = 1.16639 x 10"^ GeV-^ 

Fw = 2.06 GeV, Tz = 2.49 GeV, Vckm = 1 ■ 

The top-bottom quark-loop contribution has been evaluated using Mj = 178 GeV and 
Mb = 4.4 GeV. To study the signal-background interference we have chosen three Higgs 
mass values {Mh = 140, 170, 200 GeV), with the corresponding Higgs widths Th = 
0.008235 GeV {Mh = 140 GeV), Th = 0.3837 GeV {Mh = 170 GeV), and F^ = 1.426 GeV 
{Mh = 200 GeV), as calculated by HDECAY [48]. The weak mixing angle is given by 
Cw = M\y/Mz, = 1 — c^, and the electromagnetic coupling has been defined in the 
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scheme as aa^ = \/2G^M'y^s1^/tt. The masses of external fermions have been neglected. 
The pp cross sections have been calculated at ^/s = 14 TeV employing the LHAPDF [39] 
implementation of the CTEQ6L1 and CTEQ6M [49] parton distribution functions at tree- 
and loop-level, corresponding to Ag^ = 165 MeV and = 226 MeV with 1- and 2-loop 
running for Q;s(/x), respectively.^ The renormalization and factorization scales are set to 
My/. Fixed- width Breit-Wigner propagators are used for unstable gauge bosons. 

In Table 1 we present the total cross section and the cross section for two sets of ex- 
perimental cuts: "standard cuts" {(istd ) and Higgs search cuts {(Tbkg) as defined above. 
The results for the gluon-fusion cross section in Table 1 include the contribution from the 
massive top-bottom loop and supersede our previous calculation [20], which was based on 
intermediate light quarks only. For reference, we also show the LO and NLO quark scatter- 
ing cross sections, which have been computed with MCFM [7]. As already demonstrated 
in Ref. [20] the gg process only yields a 5% correction to the total WW cross section cal- 
culated from quark scattering at NLO QCD. When realistic Higgs search selection cuts 
are applied the correction increases to 30%. For a discussion of the renormalization and 
factorization scale uncertainties we refer to Ref. [20]. 





a{pp W*W* eui'u') [fb], LHC 


99 


qq 


""NLO 




LO NLO 




TNLO 


(^tot 


60.00(1) 


875.8(1) 1373(1) 


1.57 


1.04 


<^std 


29.798(6) 


270.5(1) 491.8(1) 


1.82 


1.06 


Cbkg 


1.4153(3) 


4.583(2) 4.79(3) 


1.05 


1.30 



Table 1: Cross sections for the gluon and quark scattering contributions to pp — » W*W* tDl'v' 
at the LHC {^/s = 14 TeV) without selection cuts (tot), with standard LHC cuts {std: pT,e > 20 
GeV, \r]e\ < 2.5, '^j, > 25 GeV) and Higgs search selection cuts (bkg) applied. The Higgs signal has 

not been included. The integration error is given in brackets. We also show the ratio of the NLO 
to LO quark scattering cross sections and the ratio of the combined NLO+gg contribution to the 
NLO cross section. 

The importance of the top-bottom loop contribution can be inferred from Table 2, 
where we compare the results based on intermediate light quarks of the first two generations 
only [20] to the contribution of the top-bottom loop and the interference between massless 
and massive quark loops. We find that the top-bottom loop increases the theoretical 
prediction by 12% and 15% for the inclusive cross section, atot, and the cross section 
with standard cuts, Cstdi respectively. After imposing Higgs search cuts, however, the 
contribution of the massive quark loop is reduced to 2% only, which is almost entirely due 
to interference with the massless loop amplitude. The reduction can largely be attributed 

®Wc observe a relative deviation of 0(10"'*) when comparing cross sections obtained with the LHAPDF 
implementation of CTEQ6 to those obtained with the original CTEQ6 implementation. 



-10- 



to the cut on A(f)j',u as can be seen in Fig. 2: while the impact of the top-bottom loop is 
sizeable in most of phase space, it is strongly reduced in the region AcpT/e < 45° selected 
by the Higgs search cuts. In Fig. 2, we also see that off-shell effects slightly decrease the 
cross section and become negligible for almost back-to-back charged leptons. Allowing for 
arbitrary invariant masses of the W bosons changes the complete background cross section 
with standard cuts by -2.8%, which increases to -6.4% when Higgs search cuts are applied. 





a{pp W*W* m'v') [fb], LHC 


quark loop 
generations 1,2 


quark loop 
generation 3 


gen. 1,2,3 
gens. 1, 2 


interference 

gens. 1,2,3 
[gens. 1,2]+ [gen. 3] 


(^totigg] 


53.64(1) 


2.859(3) 


1.12 


1.06 


cTbkgigg] 


1.3837(3) 


0.00377(2) 


1.02 


1.02 



Table 2: Cross sections for the gluon scattering contribution to pp W*W* Ivl'v' at the 
LHC (yi = 14 TeV) without selection cuts (tot) and Higgs search selection cuts (hkg) applied. We 
show the cross section with 2 massless generations [20], the contribution of the top-bottom loop, 
and the size of the interference effects. The Higgs signal has not been included. The integration 
error is given in brackets. 

An interesting distribution, which we did not display in Ref. [20] is the transverse 
mass distribution dcr/dMr, where one uses a measurable proxy for the Higgs transverse 
mass defined by Mr = {Et,u + ^tY - {pt,U + PtY, with Et,u = \JPT,u + ^li ^"^^ 

~ ''^u [^^]- Cuts on the transverse mass Mt provide an additional handle 

to suppress the background with respect to the Higgs signal, see e.g. Refs. [15,45,46]. 
In Fig. 3 wc compare the Mj^-distribution of the qq background in LO and NLO with 
the contribution from gluon-gluon scattering before and after applying Higgs search cuts. 
The figures reveal that the ghion-induced contribution becomes the dominant higher-order 
correction to the background process after Higgs selection cuts have been imposed. 

We now turn to the discussion of the interference effects between the gluon-gluon 
induced signal and background processes. Table 3 shows cross sections for the signal and 
gluon-fusion background with and without interference. We show results for Mfj = 140, 
170 and 200 GeV, spanning the Higgs mass range for which the H WW decay mode 
is of particular relevance. It turns out that the interference effects are quite small and 
never exceed 10% of the gluon-induced signal plus background cross section. Adding the 
NLO background contribution from quark scattering, the overall effect of the interference 
term is always less than 5%. These results are consistent with the small signal-background 
interference observed in pp ^ ^ 77 [51,52]. Due to the imposed jet veto, we expect only 
small effects at the LHC when NLO corrections are taken into account for the signal [53] 
and gg background. 

In Fig. 4 we finally show the W~ invariant-mass distribution without applying selec- 
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alfjgi^ 11) — W*W* — (uf'i/] [n-)] 


cut selection 


tot 


hkg 


MH[Kj(e\\ 


140 


170 


200 


140 


170 


200 


cr [signal] 


79.83(2) 


116.23(3) 


75.40(2) 


1.8852(5) 


12.974(2) 


1.6663(7) 


(7[signal + bkg(gg)] 


132.50(5) 


174.58(9) 


134.46(5) 


3.174(2) 


15.287(6) 


3.413(2) 


o-[sig+bkg{gg)] 
<r [signal] +<T [bkg (gg) J 


0.948 


0.991 


0.993 


0.962 


1.062 


1.108 



Table 3: Interference effects between the signal and gluon-induced background processes. Details 
as in Table 1. 

tion cuts for Mh = 140 GeV < 2Mw- Although this distribution is not measurable it is 
instructive in understanding interference and off-shell effects. We compare the signal and 
the gluon-induced background with and without interference effects. In addition to a pro- 
nounced, narrow resonance at Mw, the W~ invariant mass distribution exhibits a second 
small and broad resonance at about 55 GeV, which stems from the kinematic constraint 
imposed by the dominant Higgs resonance in the signal process. We note that without 
any selection cuts and for Higgs masses around 140 GeV, the gluon-induced background 
exceeds the Higgs-boson signal in the resonant region near M\y- 

4. Conclusions 

We have presented the first complete calculation of the loop-induced gluon-fusion process 
gg — — ivl'i/', including intermediate light and heavy quarks, and studied its 

importance for Higgs boson searches in the H WW channel. We find that the top- 
bottom loop, which had been neglected in Ref. [20] contributes at a level of about 10-15% 
to the inclusive gluon-induced cross section but is strongly suppressed after Higgs search 
cuts have been imposed. We have also studied interference effects between signal and 
background processes and found them to be small (about 5% or less) in the relevant Higgs 
mass range between Mh = 140 — 200 GeV. We provide the GG2WW package [40], a public 
parton- level Monte Carlo program and event generator for the process gg W*W* 
dul'v' that can be used to calculate cross sections with any set of cuts or any kind of 
differential distribution, or to generate weighted or unweighted events for experimental 
analyses. 
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Figure 2: Distribution in the transverse-plane opening angle of the charged leptons ^(^it^u for 
•pp W*W* ivi'v' at the LHC. Displayed is the background process from gluon scattering 
without (dashed) and with the top-bottom loop (solid) - the latter also without off-shell W effects 
(dotted) - each multiplied with a factor 15; the background from quark scattering at NLO (dot- 
dashed); and the signal process for Mh = 170 GeV (long-dashed) multiplied with a factor 5. Input 
parameters as defined in the main text. Standard LHC cuts have been applied. The Higgs search 
cuts select the region /S.(f)T,u < 45° left of the vertical dashed line. 




Figure 3: Distribution in the transverse mass Mr (as defined in the text) with standard cuts 
(left figure) and Higgs search cuts (right figure). Displayed are the total background from quark 
scattering at NLO and gluon- fusion (solid), and from quark scattering alone at LO (dotted) and 
NLO (dot-dashed). 
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Figure 4: Distribution in the W~ invariant mass for Mh = 140 GeV. The resonant region is 
magnified on the right hand side. Complete gg background without signal (dashed, blue), signal 
only (dotted), signal and background with interference (solid) and signal and background without 
interference (dot-dashed) are displayed. No selection cuts are applied (tot). 
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A. Auxilliary vector relations 



t-S3 U- S3 

PS = Pi -I P2 + P3 

S S 

U — Si t — S4 
P4 = Pi -\ P2 - P3 



Pi ■ P3 = P2 ■ P3 = 
K = e(Pl>P2,P3,/^) 

_2 det G 

To define the basis used in Eq. (2.6) we exploit that 

1 — S3 u — S3 

P3- J3 = Pi ■ J3 H P2 • J3 + P3 • -/s = , 

S S 

U-S4: t-S4 , 
P4 . J4 = . J4 H p2 • J4 - P3 • J4 = , 

s s 

2 s 

Js- J4 = -ipi ■ J3P2- J4+P2- J3 Pi ■ Ji) P3 ■ J3P3 ■ Ji 

S Ut — S3S4 

4 

+ -r-L ^(iPl,P2,P3, J3) ^{Pl,P2,P3, J4) , 

S{ut — S3S4) 

eipi,P2,J3,J4) = : \p3 ■ J3ko ■ Ji-P3- J^ko ■ J3] ■ 

ut — S3S4 

B. Dimension splitting formulae 

When using dimensional regularization in combination with genuinely 4-dimensional ob- 
jects, one is forced to apply a calculational scheme to deal with the 75 problem. We apply 
standard dimension splitting rules [28] and use an n-dimensional loop momentum k and 
gamma matrices ^y^, but work with 4-dimensional external momenta. The following rules 
are sufficient to evaluate the diagrams in Fig. 1: 

k — k ~\~ k J k — k ~\~ k f 

7 = 7+7 , {r,r} = o , {r,r} = 5''" , 

P^jli^ =Pjlii , Pj ■k=pj-k , 

{7^75} = , [7^75] = 0. 

All hat objects are defined in 4 dimensions, whereas the ones with tildes are (n — 4)- 
dimensional remnants. Remaining integrals that contain remnants of the n-dimensional 
algebra, i.e. factors of {k • k)" are evaluated with the following relations [54]: 



d"fc jk ■ ~kT _ ^ ^^„ r(a-e)n-4 ^ 



Tn+2a 



J iTrn/2(k^-M^)^ ^ r(l-e) 2 

r d^^k {k-kTk^k^ ^ r(a-e) n-4n + 2a 2+2a 

J iW2(fc2-^^)^ r(l-e)^ 4 n ^ 
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The n-dimensional iV-point integral is defined by [27] 

d'^k 1 



N 



iW2nf=i((fe-r,-)2-m2) 

N 



(-l)^r(A^ - n/2) / dxi...dxjv<5(l-^xO(M2)»/2-^, 



1=1 

with Tj = pi + • • • + and 



1 ^ 

= - ^ Xi^ij-Xj , Sij = -{n- rjf + ml + . 
Below we list all required integrals, expanded to the relevant order in e. 

4,a=i 7 iW2nti((A:-r,)2-m2) 8+ 



/-« = r k^ = 
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